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I I Abstract 

Ch ■ 

j^ . We present a formulation of relativistic classical and quantum me- 

^ I chanics in the Hilbert space with a linkage through the Ehrenfest 

quantization. We start with the covariant form of the Lorentz force 
and consistently derive the Dirac equation along its classical counter- 
part as a classical spinorial equation. This equation fills a missing 
Q^ I link between relativistic quantum and classical mechanics and takes 

^^ ■ the role of the Koopmann-von Neumann equation in the relativistic 

jy-s . regime. 

^ : 1 introduction 
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One of the first attempts to describe classical mechanics in terms of 
the Hilbert space was done by Koopman and von Neumann [ll |2]. 
This approach was further explored by Deoto and Mauro O [U [5] , and 
more recently enriched with the help of the Ehrenfest theorem in [6] . 
A different approach based on the formulation of classical/quantum 
mechanics in terms of the Algebra of Physical Space (APS), as a ge- 
ometric algebra variety, provided new insights on the appearance of 
the spin in a classical relativistic framework [71 [El [U [T0[ lllj . This 
work carries further the Ehrenfest quantization method described in 
[6] by developing a covariant relativistic formulation of quantum and 
classical mechanics. Our formulation starts from the covariant form 
of the Lorentz force, which allows us to deduce the Dirac equation as 
well as the corresponding equation that appears in the classical limit. 



This classical equation, which constitutes the most important novelty 
introduced in this paper, is a relativistic classical spinor equation in 
the time-extended phase-space that shows the remarkable presence of 
the interaction of the spin with the external electromagnetic field and 
a spin-orbit interaction. 

2 Background: The Lorentz force 

The Lagrangian of the relativistic particle in an electromagnetic field 

is 

^ "m ,, .,, ITT- 

C = —u''u^ + eA^'Uf, + —, (1) 



where u'^ = -^ is the four velocity also known as proper velocity and 
the Minkowski metric is taken with diagonal elements {1,-1,— 1,-1}. 
More precisely, this is the time-extended form of the Lagrangian of a 
particle with electromagnetic interaction. In this formalism the con- 
dition 

n^n^ = 1 (2) 

is not enforced as a constraint but recovered as an integration condi- 
tion in the moment to find physical solutions. 

The canonical momentum is denoted with covariant indices as 

^M = ^ = "^^/^ + e^M (3) 

Let us remember that the physical momentum and four-vector po- 
tential are denoted with contravariant indices p^ = {E, p) and A^^ = 
((/), A), where bold symbols stand for vectors in the standard Euclidean 
space. 

The time-extended Hamiltonian is 

n^P,n^-C=^P"-''^J^^-'^^^-^. (4) 

This Hamiltonian does not depend on the parameter s explicitly, so it 
is a dynamic invariant integral. The integration that leads to physical 
solutions occurs when 

H = 0, (5) 

which implies ([2]). 



(8) 



The equation of motion can be obtained accordingly with the ap- 
pHcation of the time-extended Poisson brackets defined as 

{F,G} = ——-—— (6) 

such that 

f = {F,«) (7) 

The four-velocity is recovered as 

(h^ _ dn 

ds dpfj, 

dx^ p^ — eA'^ 

ds m 

and the canonical four-force equation is 

% = e(8„A„)tt'' (9) 

ds 

The left side of this equation can be expressed as 

dpa dur, dAn dua dAn dx'^ duu dAu „ , , 

J^ = m—^ + e— ^ = m— ^ + e—^—- = m— ^ + e—^u"" (10) 
ds ds ds ds ox'^ ds ds ox^ 

and introduced in Q to obtain the covariant form of the Lorentz force 
nil 

m-^ = e{d^A, - d,A^)u'' 
ds 

m^ = eF,,u\ (11) 

where the Faraday electromagnetic tensor is defined as -Fui/ = daA 



flU '^fl-^U 



duA^. 



3 The Derivation of the Dirac Equa- 
tion 



According to ([9]) and ([8]), the dynamical equations in consideration 
are 

dpf, 

ds m' 

dx^" p^" - eAf" , , 

-r = • 13 

ds m 



^ = -id,A,W-eAn, (12) 



\\P,M - 


-- ^{d,A,W- 


-eA" 


J[^^G] - 


m 





The average of these equations over an assemble of particles written 
in the Hilbert space formalism (see, e.g., |6J) takes the form of the 
Ehrenfest theorems 

"<'''•'''''•'''''» = (*W|^(8„i„)(j.--ei^|*W>, (14) 
as m 

^<~n')\fms)) ^ (>p(,)|r^£> 

as m 

Stones 's theorem gives 

h^^^=iG\^{s)), (16) 

as 

which can be viewed as the definition of the generator of motion G. 
Following the non-relativistic case [B] , we transform the Ehrenfest the- 
orems (J12p and ()13p into the following general commutator equations 

(17) 

(18) 

Assuming that G = D{x^,p^j) and the following commutation rela- 
tions 

[x^p,] = -iM^, (19) 

we get the differential equations for the generator 

1^ = -^{d,A,W-eA^l (20) 

dD _ pf-eA^' 
dpfj. m 

The solution for the generator D has the same form as the classical 
time-extended Hamiltonian ^, which reads 

D = —{p^ - eA^iP^ - eA,) - -. (22) 

Assuming that the underlying Hilbert space is four-component, we 
can upgrade the operator D to the product of operators in a Clifford 
algebra with basis 'y^ obeying 

ji^j'^ + Y-f'' = 2g^"' X 1 (23) 



such that the generator takes the form 

D^B = ^ (^YiP^. - eA^) + rnj (^YiPu - ei.) - rnj . (24) 

This step looks more natural observing that, for commuting operators 
or classical variables, one would have D = D x 1, or in other words 
D = irr(D). 

The Dirac generator can also be written as 

D = {:^-e4 + m){^-e4-m). (25) 

2m 

The wave function equation that is consistent with the condition ([5]) 

leads to the following expression known as the quadratic form of the 

Dirac equation 

{^ - e4 + m){^ - e4 - m)^ = 0, (26) 

which is the composition of the more fundamental first order Dirac 
equation for positive and negative mass, respectively 

{^-e4- m)^ = 0, (27) 

{^-e4 + m)^ = 0. (28) 

4 The Derivation of the Relativistic 
Koopman-von Neumann Equation 

The relativistic form of the Koopman-von Neumann equation can be 
obtained starting from the general commutator equations 

i[P^,G] = ^{d^A,W-eA-), (29) 

m 

^[X^M = ^ , (30) 

supplied with the commutation relationship for the position and mo- 
mentum operators 

[X^',Pu\=^■ (31) 

Assume the dependence of the classical generator G = K{X'^, P^, A^, 0^) 
on two additional operators A^ and G'^ that satisfy the following com- 
mutation relations 

(32) 
(33) 

(34) 



[x^,K] - 


= -iS"., 


[p,,&n - 


= -^S\, 


[A„en = 


= 0. 



Then, the commutator equations reduce to the system of differential 
equations 

^ = ^(9,i.)(P^-ei^), (35) 

dK P^'-eA^' . . 

whose solution is 

iC = 1(PM _ ei'^)(A^ + ed,A^&) + /(X, P), (37) 

m 

with / as an arbitrary function of the classical operators JC^ and P^. 
Finally, the relativistic Koopman-von Neumann equation can be 
readily stated as 

Kip = 0, (38) 

with the substitution of A^, and Q'^ by the corresponding differential 
operators that preserve the commutation relations (I32p ~ (l34p . 



5 The Derivation of the Vlasov Equa- 
tion 

The classical dynamics of a relativistic particle can be expressed in 
terms of coordinates and velocities instead of coordinates and mo- 
menta, which is in fact more common, such that 

m^ = eF^.u'^ (39) 

as 

^ = n^ (40) 

as 

The application of the Ehrenfest quantization leads to the following 
commutation equations 

i[u^,G] = eF^.u-' (41) 

i[x^',G] = u''. (42) 

The rest of the derivation of the generator resembles the case treated 
in the previous section with a generator G = V{x'^,u^,A^,Q'^), such 



that 

dV 

-£ = eF,.^^ (43) 

dV 

-f = n^ (44) 

This differential equation can be readily integrated, obtaining the fol- 
lowing generator 

V = n^A^ + eF^.n'^e^ + f{x, u). (45) 

This equation can be upgraded into a different context for the descrip- 
tion of a plasma without collisions and with long range interactions. In 
this case, F^i, takes into account both the external field and the effec- 
tive field created by the particles in a self-consistent fashion. Choosing 
f{x,p) = 0, the Vlasov equation [12j for the density p in the context 
of the description of plasma reads 

u''d^p + eF^,u''^ = 0. (46) 

OUf, 

6 The Derivation of the Koopmann- 
von Neumann- Dirac Equation 

The commutator relations from the Ehrenfest quantization are 

Up^,G] = ^{d^A,W-eAn, (47) 

n m 

l[i.,G] = t^^. ,48) 

h m 

Assuming the dependence of the generator G = W{x^,p^,\^,9^) on 
the quantum operators x^ and p^ as well as on the auxiliary operators 
A^ and 0^ obeying the commutation relations 

[x'',pu] = -ihKS''^, (49) 

[x^^,X^] = -i5^^^, (50) 

[p^,n = -i6\, (51) 

[X^,r] = 0, (52) 



one reaches the fohowing differential equations 

i ( ., dW .dW\ 6/0 4 \r'~U 2v\ /roN 

nr^^-^^J = -(a,^.)(p^-e^^), (53) 

i ( ^ dW dW\ p^'-eAf , ^, 

- -ihK- 1-- = . 54 

h V op^ dXf,J m 

The solution of these equations reads 

1 ^ ^ 777 

W = ^—ip''-eA'')ip^-eA^)---fix'' + hKe'',p^-hKX^), (55) 

with / as an arbitrary function. Retracing the steps when obtaining 
the Dirac equation, we upgrade the generator W to the Clifford algebra 
spanned by the j^ matrices 

1 -^ -^ m 

W^W = -—^^^[p^-eA^)Y{Pu-eA,)---f{x^'+hKe^',p^-hK\). 

(56) 
The classical generator is obtained in the limit k — )• if / is chosen to 
be 

f = 'y^{p^-hKX^-eA^{x+hK9))^''{pu-hKXu-eAi,{x+hK9)) — -, 

(57) 
along with the following transformation of operators from the current 
quantum operators to a new set of classical operators (|34|) 

(58) 

(59) 

(60) 
(61) 

which have the remarkable property of preserving the commutation 
relations of the original quantum operators. The connection between 
quantum and classical operators are schematically represented in Fig- 
ured] 

The quantum generator, written in terms of the classical operators, 
is 

W = ^^^P^YK + —rP^.Y{At-AZ) 
m niK 



x" -- 


= x^ ^""e^ 


Pm = 




K - 
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Figure 1: Schematic layout of the quantum and classical operators. 



with i+ = A+{X^' + ^e^) and i" = 1^7(1^ - ^e'^). 



From (|62[) . we get 
leading to 



KK = hm W, 



m m 

m 2m 



(63) 

(64) 
(65) 



which can be arranged to derive the classical limit of the Dirac gen- 
erator 



K = -Y{P^. - \)Y{^^. + 5„i.e") - ^[7^A^,7M,]. (66) 
m Im 



The commutator can be manipulated as 

[7'^A^,7'^i,] = i[7^71i^/.. + [7^71^AM 



(67) 



revealing the interaction of the classical spin with the external elec- 
tromagnetic field on the left and the classical spin-orbit coupling on 
the right. 



7 General features of the Koopmann- 
von Neumann-Dirac equation 

The Koopmann-von Neumann-Dirac equation (KND) is defined in 
terms of the application of the generator (|66|) on a spinor field 

K^ = 0. (68) 

In particular, the KND equation without spin interaction can be writ- 
ten in terms of the APS representation of the geometric algebra within 
the phase space as 

K = i(P-^,(a + a„^l^). (69) 

where the spinor ^ is represented by a 2 x 2 complex matrix with the 
current and density in this representation as 

j = ^^^, (70) 

p = -rr(^'f^). (71) 

Not every solution of the KND equation represents a consistent 
classical system. A necessary condition for a truly classical solution 
is that the current needs to satisfy the KND equation 

K*^"^ = ^ ^ G classical world. (72) 

Hence, the KND equation can be reformulated as the equation for the 
current 

Kj = 0, (73) 

which constrains the form of the spinor ^. According to the singular 
value decomposition, the spinor can always be written as the product 
of a Hermitian part B and a unitary part U such that 

^ = BU. (74) 

Applying the KND operator on the corresponding current leads to 

Kj = KB\ (75) 

which implies that if a solution ^ of the KND equation is purely 
Hermitian, then the corresponding current satisfies the KND equation 
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as well. In other words, the Hermiticity of the spinor is a sufficient 
condition to be a classically valid solution 

xjrt = vi/ ^ Kj = ^ * G Classical solution. (76) 

A less restrictive sufficient condition for the current to satisfy the KND 
equation reads 

[^, ^t] = A K^t = ^ ^ G Classical solution (77) 

The KND equation must be invariant under gauge transformations 
A ^ A -\- dx- This fact can be translated into the possibility to be 
able to act on any spinor $ with the following operator 

n = e^Y>{ed,x^y (78) 

and still maintaining TZ^ as a valid solution of the KND equation. 

The KND equation is in general very difficult to solve but some 
particular analytical solutions can be found for cases involving con- 
stant electromagnetic fields. The case of a constant electric field along 
the x^ direction demands the following equation to be solve 

d + ed^^{x^)^^^ = Q, (79) 

with d^(j){x^) = —Eq. a specific spinor solution can be pursued in the 
form 

^ = f{x,p)V^, (80) 

with f{x,p) as a real scalar function and a constant Hermitian spinor 
y/u. Introducing ()80p into the KND equation, the differential equation 
for f{x,p) is 

dof -eEo^f = 0, (81) 

dp3 

which has the following solution 

f = f{po,Pi,P2,P3 + Eox^-p3), (82) 

with p3 as a constant and the shell condition for physically allowed 
solutions 

{po - Hz)f - {pif - {P2f - {psf = m' (83) 
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A particular solution in the configuration space where all the particles 
have an initial four-momentum p = po + Pi + p2 + p3 at time x^ = 
can be found by selecting 

« = -, (84) 

where the square root can be calculated with the help of formula (J112p 
such that the final solution is the following Hermitian expression 



*P = /(Po,Pi,p2,P3 + e^oa;°-p3)\/ — 

y m 



(85) 



The particle current is calculated as 



j = *# = /(po,Pi,P2,P3 + e^ox° - Pa)'-, (86) 

m 

which is also a solution of the KND equation. In this case, the density 
of particles is identified as 

p = irr(vI/M/t) = /(po,Pi,P2,P3 + eEox'> - Pa)'^, (87) 

4 m 

where po is obtained from the shell condition. If / represents a 
Dirac delta function, the solution corresponds exactly with the single- 
particle solution because the particle is constrained to be located at 

-p^ + eEox^ - Pa = 0. (88) 

For example, if the particle starts at the rest, this condition can be 
written as 

dx^ e^ .ggx 



with 7 = wl — (^)^- Further arrangement of this equation leads to 

(90) 



dx^ -0^ 



^ m 



'"° A/i + (-°f)^' 

which can be integrated to obtain a well known result 



i» 



/SJ%(^. m) 



for the initial condition x^ (0) = | ^ | • 
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8 Conclusions 

We developed a unified formulation of relativistic quantum and clas- 
sical mechanics in the Hilbert space. We applied the Ehrenfest quan- 
tization based on the classical covariant Lorentz force to obtain the 
Dirac equation as well as the Koopman-von Neumann-Dirac (KND) 
equation in the classical limit as a relativistic classical spinor equation. 
The existence of this equation has important consequences concerning 
the nature of the spin and its full consistency within the framework 
of classical mechanics. 

A Canonical commutation relations 

The commutation relation 

[x^',p,] = -id^'.h, (92) 

are satisfied with the following substitution 

xf" -^ x^ (93) 

Pf^ ^ '^^ (94) 

which is in agreement with the standard representation 

d 

E -^ ih— (95) 

dt ^ ' 

p -^ -ihV (96) 

Similarly, the following commutator relation 

[x^A,] = -^<5^, (97) 



is satisfied with 



(98) 

V - ^1;^ (99) 
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B Standard relativistic Lagrangian 

The action on both forms must the the same 

I Cds = I Ldt (100) 

/"/:^dt = f Ldt (101) 

so, up to exact differentials one has 

L = 4, (102) 

which leads to the more familiar form of the relativistic Lagrangian, 
where u^n^ = 1 is enforced as a constraint. 

C APS representation 

The APS algebra can be represented in terms of the Pauli matrices, 
such that a general expression can be faithfully represented as a 2x2 
complex matrix. The four-gradient is slitted in temporal and spatial 
components 

d = do + V (103) 

d = do-V, (104) 

with 

V = aidi+ a2d2 + a^d^ (105) 

In similar way, the four-vector potential can be represented as 

A = 9o-A (106) 

A = ao + A, (107) 

with 

A = aiA^ + a2A^ + asA^ . (108) 

For the sake of concreteness, the gradient can be fully expanded as 

9=(!'l^^ ^rf) (109) 

V5i-Fza2 do -83 J 
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The bar conjugation, also known as Clifford conjugation has the fol- 
lowing effect on any general expression 

^=f^" ^^^W^=f^}i "M (110) 

VV'21 V22J V~^21 V22 J 

The four-velocity also known as proper velocity is defined as the 
derivative respect with the proper time r 

u = — =l{l + ^), (111) 

where v is the spatial velocity. The square root of the proper velocity 

can be calculated as 

u + 1 , , 

lu = , 112 

Other important useful elements are the following complementary pro- 
jectors 

Vz = ^(1 + ^3) (113) 

Vz = ^(1-^3) (114) 

The gamma matrices in the Weyl representation can be expressed in 
terms of the Kronecker product of Pauli matrices 

7° = cJi^l (115) 

7^ = ia2 0ak (116) 

The fundamental Dirac generator can be rewritten as 

ij''{df, + eAf,)-m = ij^ do + ij'' dk + 7^ Aq + ij'' Ak - m (117) 
= CTi (g) l{ido + Aq) + Z(T2 tX" (Jkiidk + A^) (41&) 
-m td + A\ 
id + A -m J ^ ' 

The four-column spinor can be written as 

i> = {tl:^ ■ (120) 
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such that the Dirac equation is equivalent to the following coupled 
equations 

{id + A)^V3 = m^Vs (121) 

{id + A)'^^r3 = m^V^. (122) 

Conjugating the latter equation and arranging the terms we obtain 

{id + 'A)^V3 = m^Va (123) 

{-id + 1)^3 = m^^Vz (124) 

The negative sign on the second equation can be absorbed as 

id^aiP^ + A^Pa = rri^^V'i (125) 

iS^crsPs + A^Ps = m'*"^P3, (126) 

which allows to write the Dirac equation in the following compact 
form 

id^aj, +~A^ = ni^^ (127) 
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